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Plasmon spectrum and plasmon-mediated energy transfer in a multi-connected geometry 


L. Shan, E. G. Mishchenko, and M. E. Raikh 

Department of Physics and Astronomy, University of Utah, Salt Lake City, UT 84112, USA 

Surface plasmon spectrum of a metallic hyperbola can be found analytically with the separation of variables 
in elliptic coordinates. The spectrum consists of two branches: symmetric, low-frequency branch, uj < ujq/s/2, 
and antisymmetric high-frequency branch, uj > u 0 / s/2, where cc 0 is the bulk plasmon frequency. The frequency 
width of the plasmon band increases with decreasing the angle between the asymptotes of the hyperbola. For 
the simplest multi-connected geometry of two hyperbolas separated by an air spacer the plasmon spectrum 
contains two low-frequency branches and two high-frequency branches. Most remarkably, the lower of two 
low-frequency branches exists at uj -a 0, i.e., unlike a single hyperbola, it is “thresholdless.” We study how 
the complex structure of the plasmon spectrum affects the energy transfer between two emitters located on the 
surface of the same hyperbola and on the surfaces of different hyperbolas. 

PACS numbers: 73.20.Mf, 72.30.+q,78.67.n 


I. INTRODUCTION 

The two central issues of the contemporary plasmonics are 
manipulation and focusing of light on sub wavelength scales, 
and plasmon-mediated energy transfer, see recent reviews 
Refs.QUIl 

With regard to the first issue, strong confinement of opti¬ 
cal fields at small scales is accompanied by their orders-of- 
magnitude enhancement. This enhancement can be used, e.g., 
to boost nonlinear effects or to image and detecP small ob¬ 
jects. At the core of the second issue is the light-matter inter¬ 
action. The effect of plasmon-supporting interface on a group 
of emitters located nearb)P ^ 4 is twofold. Firstly, it strongly 
modifies the radiative lifetimes of individual emitters. Sec¬ 
ondly, virtual plasmon exchange facilitates the dipole-dipole 
interaction between emitters. 

A metallic strip or a wire of finite thickness can serve as 
a plasmonic waveguide. The key approach to a plasmon 
field focusing utilizes tapering, i.e. gradual narrowing of a 
waveguide towards one encPKH Field confinement at the 
end of a tapered metal wire waveguide was demonstrated 
experimentall }® 7 29 l 

Current advances in plasmonics are related to engineer¬ 
ing of progressively more complex plasmonic structures 30 35 . 
Then it is natural to extend the theoretical study of plas¬ 
monic waveguides to these complex, in particular, multi- 
connected, geometries, which contain several disconnected 
metal-air boundaries. It can be expected that these geome¬ 
tries provide additional control over plasmon fields. Studies 
of plasmons in multi-connected geometries should encompass 
calculation of the spectra of plasmonic modes as well as in¬ 
vestigation of interaction of dipole emitters mediated by such 
modes. 

Below we consider the simplest example of a multi- 
connected geometry illustrated in Fig. [T] Geometries of a 
metallic “neck,” Fig. [lk and of two tapers separated by a nar¬ 
row air “groove,” FigTjlJ), are dual to each other. Both geome¬ 
tries possess a characteristic length scale ~ a, which is the 
size of the gap or the width of the constriction. This scale, be¬ 
ing much smaller than all other scales in plasmonic structure, 
see Fig. |T]c, will therefore determine the plasmon spectrum 



FIG. 1: Simplest examples of multi-connected plasmonic structures: 
a) “neck” and b) “gap.” The plasmon spectrum of these structures can 
be found analytically if metallic surfaces are confocal hyperbolas; 2 a 
is the distance between the foci, c) Schematic view of a plasmonic 
array 36 . The field distribution in the array is determined by the points 
of contact of metallic islands. 


and the plasmon field distribution of the entire structure. 

We take advantage of the fact that the problem of calcu¬ 
lation of the plasmon spectrum in geometry shown in Fig. [T] 
can be solved exactly if the metallic surfaces are hyperbolic. 
A natural unit of momentum (wave number) is q ~ a -1 ; the 
meaning of q in the absence of translational symmetry will be 
clarified below. One advantage of the analytical approach over 
numerical methods applied to concrete sets of parameterl^is 
that an analytical solution allows to establish general proper¬ 
ties of the plasmon spectrum. 

We compare the analytical results for two hyperbolas, 
Fig. [T] to the geometry of a single hyperbola. For a single 
hyperbola, the plasmon spectrum, resembles the spectrum of 
a finite-width metallic strip. It consists of two branches: the 
low-frequency branch, uj < uj 0 /s/2, where cc 0 is the bulk plas¬ 
mon frequency, corresponds to a symmetric mode, while the 
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FIG. 2: Geometry of a single hyperbola. In elliptic coordinates £, 77 , 
defined by Eq. { 3 }, the metal with dielectric constant, e(u), occupies 
the domain (—rj 0 ,rj 0 ). Outside the hyperbola there is air with e — 1. 



FIG. 3: (Color online) Plasmon spectrum of a single hyperbola for 
two opening angles 2rj 0 = 7r/10 (red) and 2rj 0 = 7 t/3 (blue) is 
plotted from Eq. (TT) for a symmetric mode and Eq. (12) for anti¬ 
symmetric mode. At large wavenumber, m, both branches approach 
the flat surface plasmon frequency, uj 0 / \[2. 


high-frequency branch, uj > cc 0 /a/ 2, corresponds to antisym¬ 
metric mode. Adding the second hyperbola, Fig. [I] leads 
to the emergence of two additional branches in the plasmon 
spectrum. The frequency of the additional symmetric mode 
is lower than that for a single hyperbola, while the frequency 
of the additional antisymmetric mode is higher a single hy¬ 
perbola’s mode. Since the plasmons mediate the energy ex¬ 
change between the emitters close to the surface, the complex 
plasmon structure in a multi-connected geometry leads to a 
nontrivial frequency dependence of this exchange, which we 
study analytically. 

The paper is organized as follows. In Sec. II we present a 
detailed analysis of the plasmon spectrum and the field distri¬ 
bution of the plasmon modes in a single hyperbola geometry. 
This analysis is then extended to the geometry of two hyper¬ 
bolas is Sec. III. Excitation of plasmon modes by an emitter 
located at the metal boundary in a multi-connected geometry 
is studied in Sec. IV. Section V concludes the paper. 


II. PLASMON SPECTRUM OF A SINGLE HYPERBOLA 


Surface plasmons correspond to the solution of the Laplace 
equation, V (eV^) = 0, for the electrostatic potential, 

$>(x,y), which propagates along the boundary of the hyper¬ 
bola and decays away from the boundary. The variables in the 
2 D Laplace equation can be separated in elliptic coordinates 

x = a cosh £ cos 77, y = a sinh£ sin 77, ( 3 ) 


where £ is positive and 77 changes in the interval 
—7r < 77 < 7r. The Lame coefficients are the same for 

both coordinates, = \/ cosh 2 £ — cos 2 77. In the new 

coordinates the Laplace equation, 



/ <9 2 <f> <9 2 <F\ 

\ d£ 2 drf ) 


= 0 , 


( 4 ) 


has the structure similar to its usual Cartesian form. The so¬ 
lutions of Eq. ([ 4 ]) should satisfy the boundary conditions at 
77 = 770 (upper half of the hyperbola) and at 77 = —rj 0 (lower 
half of the hyperbola): the potential and the normal compo¬ 
nent of the displacement field, D n = eh^dQ/drj, should be 
continuous at this boundary. 


A. Elliptic Coordinates 


B. Plasmon Dispersion 


We start with a single-hyperbola geometry, Fig. [ 2 ] The 
boundary of a metal is defined by the equation 


/ x y 

( y ) 2 

V COS 770 / 

V sin rj 0 J 


The metal is described with a dielectric permittivity, 

e(co) = 1-^f, (2) 

C d z 

and occupies the inner part of the hyperbola, which ap¬ 
proaches the asymptotes y = db x tan rj 0 at large x. Outside 
the hyperbola there is air with e ml. 


The geometry in Fig. [ 2 ] is symmetric with respect to 
the change of sign of the 7/-axis. Correspondingly, the so¬ 
lutions of Eq. 0 can be classified into symmetric (even), 
T> s (x, y) = —y), and antisymmetric (odd), T> a (x, y) = 

— <£> a (x, —y). In hyperbolic coordinates the reversal of the 
sign of y amounts to the change 77 -A —77, so that <F S (£, 77) = 

®s{€, ~V) and $«(£, V) = -$a(£, ~v)- 
Consider first the symmetric modes. For these modes, a 
general solution of Eq. (| 4 ]) inside the hyperbola, I77 1 < 770, 
which propagates along £ and grows with 77 from 77 = 0 to¬ 
wards the boundary, has the form 


<F s (£,77) = A exp(im£) + <aexp(-ira£) cosh(77i77). ( 5 ) 
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FIG. 4 : The geometry of two hyperbolas with the opening angles 
2 rj 0 and 2771 . A dipole emitter close to the tip and polarized along the 
interface excites the dipoles located to the left and to the right from 
the tip at distances much bigger than the focal distance. 


In order to satisfy the continuity of the potential at 77 = ±770, 
the corresponding solution outside the hyperbola, \r]\ > 770, 
must have the same ^-dependence. It should also decay away 
from the boundary. This specifies the form of < 4 > s (£, 77) in the 
air: 


<^(£,77) = B exp(im^) + a exp(— im£) cosh 777(77 — I77I). 

( 6 ) 


By matching (£, 77) and ed^ s /dr] at the boundary 77 = ±770, 
we obtain the dispersion equation for the symmetric modes 


s(uj) tanh ( mr \ 0 ) coth [m( 77 — rj 0 )] = — 1 . ( 7 ) 

At first glance it appears that the constant a in Eqs. 0 - 
(6) can be arbitrary. However, from the requirement that the 
function <h s (£, 77) is finite and continuous it follows that a = 
1 . Indeed, for any a 7^ 1 the ^-component of electric field, 
= —h^d&s/d^, diverges when £ -A 0 and 77 —^ 0 , since 
the Lame coefficient vanishes there. We thus conclude that 
$*(£) cos (ra£). 

Consideration for antisymmetric modes proceeds along 
the same lines. It is not difficult to see that the solution 
cos (ra£) sinh(77777) must be discarded as discontinuous across 
the region of the x-axis, where —a<x<a and £ = 0 : the 
odd character of sinh (77777) with respect to 77 —^ —77 implies 
that it can only by multiplied by a function of £ that vanishes 
at £ = 0 . This suggests the following form of the potential, 
<F a (£, 77), for antisymmetric modes inside the hyperbola 


= Csin(m^)sinh(m7?). (8) 

Correspondingly, the potential outside the hyperbola, which is 
odd with respect to y and matches the £-dependence of Eq. 0 
has the form 




FIG. 5 : (Color online) (a) Distribution of E x (for a symmetric mode) 
and E y (for antisymmetric mode), the components of electric field 
along the x-axis, is plotted from Eqs. (T8)-(2T) for the wavenumber 
777 — 2 and the opening angle 2 rj 0 = n/ 3 . The inset shows the large- 
x oscillating tail of E x . (b) The same as (a) for two hyperbolas. 


The resulting dispersion equation for antisymmetric modes is 


e(uj) coth (777770) tanh [777(77 — 770)] = — 1 . ( 10 ) 


From Eqs. 0 and ( fT()| ) we can express the plasmon frequen¬ 
cies in terms of the dimensionless wavenumber, 777, 


{ \ 1/2 

Sinh(77777 0 ) COSh[777 ( 7T - 77o)] I 

-\- } > (U) 

smh(7777r) I 


{ \ 1/2 
cosh(m7?o) sinh[m(7r - T] 0 )] I 

-- 1 ' m 


Examples of dispersions Eqs. ( [IT] ), ( fl2| ) are shown in Fig. [ 3 ] 
Below we list some general properties of these dispersions: 

(i) A hyperbola reduces to a plane for 770 = tt/ 2 . Then 
Eqs. (TT) , (T2) reproduce the expected result, uj s (m ) = 

(m) = cCq/V^, for the frequency of a dispersionless sur¬ 
face plasmon. 

(ii) The spectra Eqs. (TT)-(T2) of the symmetric and anti¬ 
symmetric modes satisfy the “sum rule” relation, 


( 9 ) 


W^( to ) +w 2( to ) = uj 2 _ 


$,*(£,77) = Dsin(m£) sinhm( 7 r — \r)\). 


( 13 ) 
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(iii) In the long-wavelength limit, m —>> 0 , the threshold 
frequencies are 


.( 0 ) 




(14) 


The fact that cc s ( 0 ) goes to zero for small r] 0 could be ex¬ 
pected, since at small rj 0 the hyperbola is effectively a metal¬ 
lic layer with zero thickness, for which the dispersion of the 
longitudinal surface plasmon does not have a threshold 38 . On 
the other hand, uj a approaching uj 0 for r ) 0 —» 0 can be inter¬ 
preted by noticing that the oscillations of the electron density, 
accompanying this plasmon, are normal to the surface; such 
oscillations must have the bulk-plasmon frequency ccq • 


(iiii) Another notable property of the dispersions Eqs. ( 11 )- 
© is their duality , namely, 



w s (7r - r? 0 ,m) = oj a (rio,m). (15) 

Qualitative interpretation of this relation can be given for 
small r] 0 , when Eq. © relates the spectrum in of a sharp 
metal “edge” with the spectrum of a narrow “funnel”: it sug¬ 
gests that the frequency of a soft symmetric plasmon of the 
“edge” coincides with the frequency of an antisymmetric plas¬ 
mon in a “funnel.” Indeed, once the plasmon frequency is low, 
the absolute value of e(uj) is big. This, in turn, implies that the 
normal component of the oscillating electric field in the metal 
is small. In the funnel geometry, this small field amplitude can 
be realized only when opposite charges accumulate on the two 
surfaces of the funnel. The reason is that, for small-angle fun¬ 
nel, such antisymmetric arrangement is similar to the charges 
of a parallel-plate capacitor; this distribution of charges en¬ 
sures that the field of the parallel-plate capacitor does not ex¬ 
tend outwards. On the other hand, as we mentioned above, 
the low-frequency plasmon in a sharp-edge geometry is lon¬ 
gitudinal, which corresponds to symmetric amplitudes of the 
charge-density fluctuations for the two surfaces. 


C. Field Distribution 

In addition to the spectrum, cc(ra), it is instructive to look 
at the spatial distribution of electric field of the two plasmon 
modes. For a given frequency, u, the value of m is found by 
equating it to either u s (m ) 9 when c j < uj 0 /y/ 2 , or to c o a (m), 
when uj > cc 0 /\/ 2 . The found value of m(cj) is then substi¬ 
tuted into the potential distribution 4 > s — cos(ra£) cosh(ra?7) 
or <F a = sin(ra£) sinh (mrj), from which the electric field is 
subsequently calculated. 

To clarify the physical meaning of the parameter m we 
rewrite the potential 4 > s at large distances from the coordinate 
origin in the Cartesian coordinates 


FIG. 6: (Color online) Density plot of the field intensity of the plas¬ 
mon modes for a single hyperbola with the opening angle 2rj 0 = 
7 t/3 . The upper row corresponds to a symmetric mode, for which 
the potential distribution is determined by Eqs. ^ and {6|. The 
lower row corresponds to an antisymmetric mode with the potential 
described by Eqs. {8} and Left, central, and right panels cor¬ 
respond to the wavenumbers, m — 0.1, m — 0.4, and m — 0.7, 
respectively. With increasing m the plasmon frequencies approach 
cjo/a/ 2, while the field concentrates near the metal-air surface. 

The ratio of the amplitudes of 4 > s at the boundary, y = 
x tan r] 0 , and along the x-axis is equal to cosh ( mp Q ). 

The potential distribution Eq. (I 6 | ) applies inside the metal. 
In the air, this distribution differs from Eq. ( [T 6 | ) by the 

replacement arctan^^ by tt — arctan^^ and C by 

C cosh(m?7 0 )/cosh [rrt( tt — rj 0 )]. 

Consider now a symmetric plasmon propagating with a 
wavevector, q , in a metallic film of a constant thickness, 2 d. 
For this plasmon the ratio of the potentials at the boundary 
and at the center is equal to cosh (qd). Identifying d with the 
length of the arc between the x-axis and the boundary, pp 0 , 
where p = ^ x 2 + y 2 , allows one to specify the local value of 
the wavevector 

q(p) = -• ( 17 ) 

P 

With the p-dependent wavevector given by Eq. fTT] ), one 
would expect for the plasmon phase a “semiclassica” value 
dp'q(p') = mln(p/a), which is indeed the case, as fol¬ 
lows from Eq. ( [T 6 ] ). 

For a symmetric mode, the electric field on the x-axis is di¬ 
rected along x. The domain x > a on the x-axis corresponds 
to r] = 0 , so that T> s (£) = C cos (ra£). In this domain the 
behavior of electric field with x = a cosh £ has the form 


®s(x,y)\x,y»a = C COS 

x cosh 


In 


2 yjx 2 + y 2 


m arctan — 
x\ 


( 16 ) 


E x (x > a) = Cm 


sm 




\Jx 2 — i 


( 18 ) 


The field is finite at x = a, and falls off as 1/x for x a. 
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FIG. 7: (Color online) Density plot of the field intensity of the plas- 
mon modes in the geometry of two symmetric hyperbolas with the 
same opening angle, 2r] 0 = 7t/3, as in Fig. [6] The upper row cor¬ 
responds to a symmetric mode for which the potential distribution 
is determined by Eqs. ( [22} and {27}. Left, central, and right panels 
correspond, respectively, to the wavenumbers, m = 0.1, m = 0.4, 
and m = 0.7, the same as in Fig. [6] While the fields of individual 
hyperbolas are disconnected along the x-axis, they overlap along the 
y- axis. The lower row corresponds to antisymmetric modes with the 
same 777-values. The field of individual hyperbolas overlap, predom¬ 
inantly, along the y- axis. 


In the domain a cos r ] 0 < x < a we have £ = 0 and <F S = 
C cosh (77177). Differentiating with respect to x = a cos 77, we 
obtain 




FIG. 8: (Color online) Upper panel: Comparison of the plasmon 
spectra in the geometry of two hyperbolas for two values of the open¬ 
ing angle. For rj 0 = 7 t /6 < 7 r /4 the frequencies of both symmet¬ 
ric modes (blue) are smaller than the flat surface plasmon frequency 
uj Q / \/2, while the frequencies of both antisymmetric modes (red) are 
bigger than uj 0 / y/2. The relative signs of the oscillating charge den¬ 
sity along the metal surfaces are schematically illustrated to the left 
of the graph. For 770 = 27t/5 > 7 t /4 (lower panel) the positions 
of the upper symmetric mode and lower antisymmetric mode with 
respect to ce 0 invert. 


E x (a cos 770 < x < a) = Cm 


sinh ^777 arccos ^ 
yj a 2 — x 2 


( 19 ) 


Finally, in the air, in the domain 0 < x < acosr ] 0 the x- 
component of the field is given by 


= —mC I 


E y (0 < x < acosr ] 0 ) 
sinh(m?7 0 ) \ sinh [m( 7 r - arccos ^)] 


sinh [777 ( 7 r — rj 0 )] 


V a 2 — x 2 


( 21 ) 


= —mC 


cosh(m 77 0 ) 
cosh [m(ir — r ] 0 )] 


E x (0 < x < acos 77 0 ) 
sinh [m( 7 r — arccos ^)] 


Vo 


( 20 ) 


The behavior of E x (x) and E y {x) is illustrated in Fig. [ 5 ^. It 
is seen that the oscillating tail emerges only at large distance 
20 a from the tip. 


In Eqs. ( fl8| ), ( fl 9 | ) the frequency and the 770-dependence of 
electric field is incorporated in m (cc, 770), defined by the con¬ 
dition cc s (m) = u), where u s (m) is given by Eq. ( 11 ). 


For an antisymmetric mode, the electric field on the x- 
axis is directed along y. Inside the metal, where < 4 > a = 
C sin(ra£) sinh(77777), th e ^-dependence of E y is the same as 
the ^-dependence, Eqs. ( [ 20 } , of E x for the symmetric 
mode. In the air, the x-dependence of E y is different from 
Eq. ( [ 20 } and has the form 


III. GEOMETRY OF TWO HYPERBOLAS 
A. Splitting of the Plasmon Spectrum 

Consider now the geometry of two hyperbolas, Fig. | 4 } The 
boundary of the first hyperbola is defined by the same Eq. 0; 
in addition, metal occupies the domain < 77 < 7 t. The £- 
dependence of the potential for the symmetric modes remains 
the same, cos The 77-dependence, which decays away 
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from both boundaries, is characterized by four constants, 





A\ cosh(m?7), 0 < r] < r ] 0 

B\ cosh(m?7) + B 2 cosh 771(71 — 77), r ] 0 < r] < tt — r ] 1 
A2 COShm( 7 T — 77), 7 T — 77 ! < 77 < 7 T, 


( 22 ) 


These constants, Ai, A 2 , B\ and B 2 are related by continuity 
of $> s (rj) and ed^ s /dr] at 77 = r ] 0 and 77 = tt — rj 1 . These 
continuity conditions read 


-A2 cosh( 7 n? 7 1 ) = B\ cosh 777(71 — 77J + .82 cosh( 77777 !), 

—e(lS)A 2 sinh( 77 i 77 1 ) = sinh 777(71 — 77!) — B 2 sinh( 7 ?i 77 1 ). 

( 24 ) 


A\ cosh( 77177 o) = B\ COSh( 777 / 77 o) + B 2 COSh 777/(71 — 770), 
e(uj)Ai sinh(77i77 0 ) = B\ sinh(77777 0 ) — B 2 sinh 777(71 — rj 0 ), 

( 23 ) 


These relations, together with the explicit form of e(u), given 
in Eq. lead to the following characteristic equation 


sinh 2 mr ] 0 sinh 2 mri 1 


fuj 2, 1\ . 1 . , /or 1\ 

( —77-smh 77771 H— smh 777(71 — 2 t 7 0 ) —-smhi 

W 2/ 2 v ,oy Vu>2 2/ 


r 


— sinhm( 7 r — 2 %) 


( 25 ) 


The two brackets in the right-hand side describe the plasmon 
dispersion Eq. (TT] ) for a symmetric mode, while the left-hand 
side describes the coupling of the two plasmon branches. The 
two plasmons decouple when rj 1 is small. Then the dispersion 
of the upper symmetric branch is simply uj s Eq. ( 11 ). 

In order to find the dispersion of the lower symmetric branch, 
it is sufficient to set uj = 0 in the first bracket. This yields 


U} s ( to ) 


771 <1 


= CUn 


mrjx tanh 777(71 — r ] 0 ) 


1/2 


( 26 ) 


The general expression for the dispersion of the two cou¬ 
pled symmetric branches reads 


uf ( to ) = 


CCn 


[2 sinh( 777 jr)] 


1/2 


sinh( 77 i 77 o) cosh771(71 — r] 0 ) + sinh^^) cosh771(71 — 77J 

1/2 


± 


sinh 2 m ( r ] 0 — 77J cosh 2 m (tt — r ] 0 — 77J + sinh ( 2 mr] 0 ) sinh (2 mr ] 1 ) 


\ j 

1 1/2 1 


( 27 ) 


It is easy to see from Eq. ( [ 27 ] ) that for large m 1 both fre¬ 
quencies approach the surface plasmon frequency cc 0 /y/ 2 , as 
in the case of a single hyperbola. The reason is that the short- 
wavelength plasmon is “local,” a metal surface is locally flat, 
and the presence of the second surface is of no consequence 
to the spectrum in this limit. The behavior of uof (m) in the 


limit of long wavelengths is remarkable 


■(m) 


: TYlUJr) 


ra—>-0 


TloVi (?r - % - »?i) 


V0+V1 


1/2 


fVo + m \ 1/2 


( 28 ) 


The fact that ( 0 ) is determined by the “net” angle (rj 0 + 77J 
is consistent with the result Eq. ( [ 14 ] ) for a single hyperbola. 
As this net angle approaches 7r, the portion of air in this limit 
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becomes small, and (0), approaches the bulk plasmon fre¬ 
quency. The acoustic behavior of uj~ (m) is related to the fact 
that, unlike for a single hyperbola, in the geometry of two hy¬ 
perbolas a low-frequency plasmon is not reflected from the 
tip, but goes “through” the gap into the second hyperbola. 

For two identical hyperbolas, p 0 = rj l9 Eq. ( [27] ) simplifies 
to 




= u. 


( smh(mp 0 ) ^cosh [m(7r — r] 0 )} =b cosh(ra77 0 )^ | 


i 


sinh(m7r) 


} 


(29) 


Similar derivation for the antisymmetric plasmons yields 


w±(m) 

coslifm//,,) ^sinh [m( n — rj,,)] ± sinh (mr/, 
sinh(m7r) 

(30) 

From Eqs. ( [22] ), ( [30} one can trace the evolution of the plas¬ 
mon spectrum with increasing rj 0 . For rj 0 <1 the frequencies 
of both symmetric plasmons are low: 




u s ( to ) 


~ cu 0 

. ( mir \ 
mp o tanh 


V 2 / 


1/2 


(31) 


"(m) 



mrj 0 

« ^0 

V o<l 

tanh (^f) 


1/2 


(32) 


while the frequencies of antisymmetric plasmons are close to 
uj 0 . As t] 0 increases and achieves the value 7 t/ 4, the branches 
cc+(m) and collapse into a single frequency uj q /V2 

and become flat. As p 0 increases further above 7 t/ 4, the 
branches invert: cu+(m) is pushed above uj 0 /V^ and uj~(m) 
drops below it. This evolution is illustrated in Fig. [8] 

For the geometry of two hyperbolas there is a duality rela¬ 
tion, 



7o, 


t) =u a (Vo ,m). 


similar to Eq. © for a single hyperbola. 


(33) 


B. Comparison of the field distributions for two geometries 

Hybridization of plasmon fields of individual hyperbolas in 
the geometry of two hyperbolas is illustrated in Fig. [5] Fig. [6] 
and Fig. [7] The curves in Fig. [5] suggest that hybridization 
along the x-axis is rather weak and becomes progressively 
weaker while the wavenumber increases as the frequency ap¬ 
proaches uj q / y/2. This behavior is natural, since, the closer 
the frequency is to that of the flat surface plasmon, the more 
localized is the plasmon field near the metal-air interface. The 
density plot of the modes of individual hyperbolas is shown in 
Fig. [6] Compared with the latter, the double-hyperbola plot of 
Fig. ^demonstrates that hybridization of individual modes has 
a “ring”-like character for symmetric modes and the “needle”- 
like character for antisymmetric modes. Such different nature 
of hybridization can be interpreted with the help of the pat¬ 
terns of oscillating surface charges in Fig. [8] For the bottom 
symmetric mode with low frequency, the positive and negative 
charges are separated by air, whereas the dielectric function of 
the metal is large. This expels the electric field lines from the 
inside of the metal and forces them to go through the air. The 
field is strong along the y- axis, where it is parallel to the x- 
axis. 

For the antisymmetric mode, the upper and lower sides of 
the metal have opposite charges, while the frequency is close 
to cc 0 , so that the dielectric constant of the metal is small. 
The force lines of electric field are localized inside the metal. 
Thus, the two metal edges can be viewed as the plates of a 
parallel-plate capacitor. Correspondingly, the concentration 
of electric field near the x axis is analogous to the fringe field 
outside a parallel-plate capacitor. 

In Sect. II we have established that the field of the plas¬ 
mon modes at large distance, p a, from the origin be¬ 
haves as 1/p. Hence, the field intensity behaves as 1/p 2 , 
i.e. it strongly diverges down to the distances ~ a. This, 
however, does not translate into a strong enhancement of the 
net energy, J drE 2 (r), which grows only logarithmically, 
oc ln(A/a); the upper cut-off is provided by the wavelength 
of light, A = 27 tc/uj, with the same frequency uj. 


C. Two Co-directed Hyperbolas 


For completeness, in this subsection we will analyze the 
plasmon spectrum in the geometry of two co-directed hyper¬ 
bolas. Assume that the metal occupies the region p 0 < r\ < 
rji, while the regions 0 < p < p 0 and p 1 < p < tt are occu¬ 
pied by air, see the inset in Fig. [9] A straightforward general¬ 
ization of Eq. ( [25] ) gives 


sinh 2mr] 0 sinh 2 m ( 7 t — %) 


— ^"jsinhraTT — ^ sinhm(7r — 2 rj Q ) ^sinhrmr + ^ sinhm(7r — 2p x ) 

\ CU n Z / Z V UJ n z / z 


■ ( 34 ) 


4 
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FIG. 9: (Color online) Illustration of the plasmon spectrum in the 
geometry of two co-directed hyperbolas. Blue and red curves corre¬ 
spond to symmetric and antisymmetric plasmon modes, respectively. 
The spectrum is plotted from Eq. ( [34] ) for rj c = Sr] — 7 t/6 . 


To analyze the plasmon dispersion, we introduce the average 
opening angle and the ’’thickness” of the tip, 


Vc = — 2 —> Sv = Vi~ Vo ■ (35) 

An example of the plasmon spectrum with the two co¬ 
directed hyperbolas is shown in Fig. [9] In this inverted geom¬ 
etry, symmetric and antisymmetric modes interchange, com¬ 
pared with Fig. [8jb). The symmetric branches in the long- 
wavelength limit are 

/5n\ 1/2 

w+(0)=w o (—J , uJ(0)=uj o . (36) 


For the upper symmetric mode, uif (0), the signs of the os- 
cillating charges on the opposite sides of each “sleeve” are op¬ 
posite and the electric field lines are confined inside the metal. 
For this, the dielectric constant of the metal must vanish. For 
the lower symmetric mode, cc<r(0), the signs of charges on 
the opposite sides of each sleeve are the same. The field lines 
mostly stay in the air and do not penetrate into the metal. This, 
on the other hand, implies that e(uj) is large, and, correspond¬ 
ingly, the frequency of the mode is small, vanishing in the 
limit of a very thin “coating,” Sr] -4 0. 

The acoustic mode, generic for multi-connected geome¬ 
tries, is now found in the antisymmetric part of the spectrum, 
(m) -4 0 as m -4 0. The second mode, 

u;-(0)=cu o (l-^J . (37) 

lies above cc 0 / y/2 and its frequency increases with decreasing 
Sr]. 



FIG. 10: (Color online) Geometric factor for the energy transfer 
is plotted as a function of frequency from Eq. ( [43] ) for the val¬ 
ues of the opening angle 2r]± = 0.2n, 0.247T, 0 . 287 T (a), and 
2 771 = 0 . 567 T, 0 . 67 T, 0.647T (b). The value of the opening angle, 
2?7 o, is chosen 2r] 0 = 47 t/ 5 . 


IV. INTERACTION OF TWO EMITTERS AT THE 
METAL-AIR SURFACE 


Numerous studies have demonstrated that, due to proxim¬ 
ity to the metal, the lifetime of the emitter can be dominated 
by the excitation of the plasmon modes 6 ®. For a metallic 
nanoparticle plasmons significantly shorten the lifetime when 
the emitter frequency is close to the frequency of the plas¬ 
mon dipole (/ = 1) oscillations. A less trivial finding 1314 
is that the coupling to plasmons can dominate the lifetime 
when the emitter is located in the proximity to a nanowire 
with a subwavelength radius in which the plasmon spectrum 
is continuous. When the emitter is positioned close to the tip, 
the plasmon-induced shortening of the lifetime is even more 
pronounced^. This is caused by the field enhancement near 
the tip. In this regard, to calculate the emitter lifetime in 
a multi-connected plasmonic structure, we can use the same 
general scheme as developed in the previous studies. 

Consider for concreteness, a dipole emitter with frequency 
u), located at the metal-air interface, and polarized parallel to 
the interface, as in Fig. [4] The structure of the plasmon spec¬ 
trum in this calculation is captured by Green’s function for 
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electric field, 


G&mt'rf-, <*>) = £ 


uj — uj{m) 


(38) 


With the help of the continuity conditions Eqs. ( [23] ), (24) this 
ratio can be cast in the form 


1 -2 


where 77 ; ra) is the normalized tangential component of 
the electric field of the plasmonic mode and the summation 
is taken over all four modes; u(m) is the plasmon spectrum 
found above, The decay rate is proportional to the imaginary 
part of the diagonal value G(£, r/*; cc) taken at the po¬ 

sition of emitter. As the summation in Eq. ( [38] ) is performed 
over acoustic as well as optical modes, a spike-like feature oc¬ 
curs in the uj- -dependence of the decay rate. The origin of this 
feature is the divergence of the density of the plasmon modes, 



dm S [uj — cc + (ra)], 


(39) 


near the threshold frequency cc+( 0 ). 

Consider for simplicity a symmetric geometry, rj 0 = r] 1 < 
7 r/ 4 . Expanding Eq. (22) at small m, we get 


ut( m ) 

ut(0) 


m—> 0 





(40) 


where w+(0) = {2t} 0 /'k) 1 ^ 2 . Substituting Eq. 
one finds 


(40) into Eq. 


9 + M = - 

7r 


(f -Vo) (f -T)o)ut(0) (w-wt(0)) 


1/2 


(41) 

The inverse square-root singularity in the density of states 
translates into a minimum in the emitter lifetime. Similar min¬ 
imum is expected near uj = ( 0 ) and near the bulk plasmon 

frequency, uj 0 . Near uj = uj 0 / a/2 the minimum should be even 
stronger, since the divergence of the plasmon modes diverges 
as \uj — cCq/V^I -1 - 

The scheme of calculation of the plasmon-mediated energy 
transfer between the emitters near the metal-air interface is 
also well established, see e.g. Ref. [l0]for two emitters near 
a nanoparticle and Ref. 39 for donor and acceptor above the 
graphene. Here we emphasize the specifics of this transfer in 
the geometry of two hyperbolas. In particular, we study how 
the transfer rate between the red emitter near the tip in Fig. [4] 
and the blue emitter located on the left hyperbola differs from 
the transfer rate between the red emitter and the blue emitter 
located on the right hyperbola, when both blue emitters have 
the same coordinate, £. 

Quantitatively, the relation of the two transfer rates is de¬ 
scribed by the ratio of the tangential components of electric 
field at the two interfaces. From Eq. ( ]22| ) we find 


T= ( E((rj 0 ,uj) \ 


B\ cosh(ra77o) + coshra( 7 r — r\ 0 ) 
B\ coshm( 7 r — 771 ) -\- B 2 cost/rar^) 


2 


(42) 


T = 


e(uj) tanh(m? 7 1 ) sinhm(? 7 Q —r/J+cosh m(rj 1 — rj 0 ) 


(43) 


where ra(cc) is determined by the dispersion equation (21). In 
Fig. [TO] we plot the factor T as function of frequency for a 
fixed opening angle of the right hyperbola, 2r] 0 = 47 r/ 5 , and 
several opening angles of the left hyperbola. Overall, we see 
that at uj = 0 the geometric factor T rapidly increases with 
r] 1 . This can be understood from the electrostatics of ideal 
metals. Indeed, at uj = 0 the field does not penetrate into the 
metal at all. Consequently, the sharper the left hyperbola is, 
the stronger the field near the left tip becomes, resulting in 
smaller values of T(0). This, in turn, means that the energy 
transfer from the red emitter in Fig. [4] to the left blue emitter 
happens faster than the transfer to the right blue emitter. 

We also see that the behavior of T(u) is different for small, 
Fig.fTOfc), and large, Fig.fTOfr), values of r] 1 . When geometrical 
factor T(0) is small, T{uS) falls off with frequency, suggesting 
that the left and right hyperbolas become effectively “discon¬ 
nected”. To the contrary, for larger starting values T(0), the 
finite-frequency T(uS) grows with frequency, suggesting that 
the energy transfer becomes more symmetric. The crossover 
from decay to growth takes place at r] 1 « 0.197r. This value 
can be related to the peculiar behavior of the velocity of the 
low-frequency plasmon Eq. ( [28] ). Namely, this velocity has a 
minimum as function of r] 1 . By setting duj~/dr] 1 = 0, we find 
the position of this minimum at, 


fji = 0™7o) 1/2 - Vo- 


(44) 


For T] 0 = 27t/ 5, Eq. ( [44] ) yields fj 1 ~ 0.237T, 

close to the crossover value of r] 1 , which corresponds to 

dT(uj)/duj |^o=0. 


V. CONCLUDING REMARKS 


(i) The prime qualitative finding of the present manuscript 
is that, for a singly-connected geometry, the opening angle of 
the hyperbola, 2 r] 0 , defines the minimum frequency, given by 
Eq. below which the plasmon bound to the tip of the 
hyperbola does not exist. By contrast, in a multi-connected 
geometry a symmetric plasmon mode, concentrated near the 
region of the closest contact of the surfaces, exist at arbitrary 
low frequency. 

(ii) Originally, the enhancement of the electric field of a 
plasmon-polariton as it approaches a wedge, has been demon¬ 
strated analytically in Ref. [15] Calculation in Ref. ES was 
carried out in polar coordinates, i.e. neglecting the rounding 
of the wedge near the tip. This curving has been emulated by 
introducing a cutoff length ~ a. At distances much greater 
than a the form of the plasmon field in Ref. [15] is consistent 
with Eq. ( [16] ) with one important difference that the plasmon 
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field Eq. ( [T6| ) is not traveling, as in Ref. ! 13J but is a stand¬ 
ing wave instead. This structure is enforced by the boundary 
condition at the tip. 

The fact that the plasmon mode near the tip is a standing 
wave can be viewed from the perspective of focusing of light 
on subwavelength scale, when the plasmon-polariton with TM 
polarization (magnetic field along the 2 -axis) is excited at 
large (compared to the wavelength, A) distance, p , from the 
tip 18 . In the geometry of a single hyperbola this polariton 
will be fully reflected. For small enough opening angles, 2rj 0 , 
the transformation of the polariton into small-p standing plas¬ 
mon mode can be traced analytically within the semiclassical 
(WKB) description 13 Within that description, the 
metal strip at distances p^> a from the origin is replaced with 
a planar film with a thickness d p , equal to the arc distance, 
2pp 0 = d p , between the two metal surfaces. Then the semi¬ 
classical expression for the incident and reflected fields has 
the form g” 1 / 2 exp [±i f dp q p ], where q p (oj) is the disper¬ 
sion of the symmetric waveguide mode propagating along the 
film. This dispersion law satisfies the equation 


the tip for a symmetric plasmon. This ratio can be found from 
Eq. ( p~9] >. The origin, x = 0, corresponds to p = tt/2, while 
the tip corresponds to x = a cos rj 0 . Then Eq. ( [T9] ) yields 


E x (0) sin tj 0 sinh [|m(cc)] 
E x ,ti P sinh [0r - r] 0 )m(uj)\ ' 


(47) 


As the frequency grows from the threshold value cc s (0) = 

to the surface plasmon frequency cj 0 /\/2, the 

“wavenumber,” m, changes from m = 0 to m = oo. Then 
the ratio Eq. © falls off from 7rsin?7o/[2 (w — rj 0 )} mono- 
tonically. At frequencies close to uj 0 /\/2 the ratio Eq. (47) 
approaches zero, since the plasmon field is strongly localized 
near the metal surface. Such a strong change of the field is 
revealed by the exact solution in elliptic coordinates demon¬ 
strated in the present paper. 

(iii) In this paper we assumed that the plasmon field does 
not depend on z. For a general ^-dependence, <f> oc exp (m 2 ), 
the Laplace equation takes the form 


5 tanm 





The relevant question is, how close to the threshold frequency 

lj 0 (rjo/n'j is the semiclassical description applicable, i.e. 

when does the condition dq p /dp <C q 2 apply? At small fre¬ 
quencies we can replace e(cj) by — gc 2 /gc 2 and simplify Eq. 

@ ) using the identification q p — m/p , Eq. (17). Then Eq. 
) reduces to the following equation for m 


tanh 




*oVy / 2 

c 2 m 2 ) 


„V\V 2 




V c z m z / 


(46) 


and the semiclassical condition is satisfied provided that m >> 
1. As the distance p increases, the solution of Eq. grows 
starting from m = , which is consistent with Eq. Jilt. 

Thus, for semiclassical description to apply at all p bigger than 

a, it is necessary that this minimum m exceeds 1 , i.e. uj 

1 /2 

uj 0 i 70 . The latter condition suggests that the frequency of the 

incident light, being much smaller than cc 0 , should still not be 
very close to the threshold frequency. 

Another approach to the “delivery” of the light energy to 
the tip was proposed in Ref. [21] Namely, one can coat a con¬ 
ical tip of a glass fiber by a silver layer, the geometry similar 
to the one shown in the inset in Fig. [9] The wavelength of a 
plasmon in a silver layer increases away from the tip, and at 
a certain distance matches the wavelength of the waveguided 
light propagating towards the tip. At this point the light trans¬ 
forms into the plasmon and heads towards the tip. Experimen¬ 
tal realization of a coupling based on this idea was reported in 
a number of papers, see Ref. 37 and references therein. 

(ii) Naturally, the structure of the plasmon field at small dis¬ 
tances from the tip is not captured in polar coordinates. Mean¬ 
while, this structure is quite nontrivial. To illustrate this, con¬ 
sider the ratio E x (0)/E Xj ti P of the fields at the origin and at 


2 + ft 2 & 2 (cosh 2 f - COS 2 7 ])$ = 0. (48) 

Separation of variables, leads to the following equation for the 
potential’s ^-dependence, 

d 2 $ 

- 7 — r - — (m 2 + k 2 cl 2 cos 2 = 0. (49) 

OT] z 

Our results apply for the wavenumbers that exceed some 
minimum value, m m min . This minimum value can be es¬ 
timated from the observation that k can not be less than H~ 1 , 
where H is the “height” of the hyperbola, see Fig. [l] This 
determines, m rniri ^fta^a/77<Cl. 

(iiii) A finite scattering rate of conduction electrons, 7 , sets 
the maximum distance from the origin, p max , where our re¬ 
sults apply. The solutions obtained and discussed above decay 
beyond p rnax . This distance can be estimated by noting that a 
finite 7 is taken into account by replacing uj 2 uj (lj + i/ 7 ) in 
the dielectric function, Eq. 0- As a result, it leads to a finite 
imaginary part of m: Imm « 7 (dc j/9m) _1 , where duo/dm 
is the slope of the plasmon dispersion. The plasmons attenu¬ 
ate at distances where (Imm ^ 1. Since at large distances £ 
depends logarithmically on p, see Eq. ( [T6| ), we conclude that 
the maximum distance is p max ~ a exp [1 /Im m \. 
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